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THE DENSITY OF PRIMES DIVIDING A TERM IN THE 
SOMOS-5 SEQUENCE 

BRYANT DAVIS, REBECCA KOTSONIS, AND JEREMY ROUSE 


Abstract. The Somos-5 sequence is defined by ao = ai = <12 = <13 = <14 = 1 
and a-m = am - iam - 4 + a m- 2 am .-3 f or m > 5 We relate the arithmetic of 

a m—5 ~ 

the Somos-5 sequence to the elliptic curve E : y 2 + xy = x 3 + x 2 — 2x and 
use properties of Galois representations attached to E to prove the density of 
primes p dividing some term in the Somos-5 sequence is equal to //A, . 


1. Introduction and Statement of Results 

There are many results in number theory that relate to a determination of the 
primes dividing some particular sequence. For example, it well-known that if p is 
a prime number, then p divides some term of the Fibonacci sequence, defined by 
Fq = 0, Fi = 1, and F n = F n _ 1 + F n _2 for n > 2. Students in elementary number 
theory learn that a prime p divides a number of the form n 2 + 1 if and only if p = 2 
or p = 1 (mod 4). 

In 1966, Hasse proved in |Jj that if 7r even (x) is the number of primes p < x so 
that p\2 n + 1 for some n, then 

,. Teven(^) 17 
hm —(— = —. 

m —>00 TT[X) 24 

Note that a prime number p divides 2 n + 1 if and only if 2 has even order in . 

A related result is the following. The Lucas numbers are defined by Lq = 2, 
Li = 1 and L n = L n _ 1 + L n _ 2 for n > 2. In 1985, Lagarias proved (see j9j( and 
[lOl l that the density of primes dividing some Lucas number is 2/3. Given a prime 
number p , let Z{p) be the smallest integer m so that p\F m . A prime p divides L n 
for some n if and only if Z(p) is even. In [2], Paul Cubre and the third author 
prove a conjecture of Bruckman and Anderson on the density of primes p for which 
m\Z(p), for an arbitrary positive integer m. 

In the early 1980s, Michael Somos discovered integer-valued non-linear recur¬ 
rence sequences. The Somos-fc sequence is defined by Co = C\ = • • ■ = Ck-i = 1 
and 

+ Cm—2Cm—(fc—2) + ' ' ' + c m - J C m 


for m > k. Despite the fact that division is involved in the definition of the Somos 
sequences, the values c m are integral for 4 < k < 7. Fomin and Zelevinsky [3] show 
that the introduction of parameters into the recurrence results in the c m being 
Laurent polynomials in those parameters. Also, Speyer m gave a combinatorial 
interpretation of the Somos sequences in terms of the number of perfect matchings 
in a family of graphs. 
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Somos-4 and Somos-5 type sequences are also connected with the arithmetic 
of elliptic curves (a connection made quite explicit by A. N. W. Hone in [5], and 
E3)- If a n is the nth term in the Somos-4 sequence, E : y 2 + y = x 3 — x and 
P = (0,0) £ -E’(Q), then the ^-coordinate of the denominator of (2n — 3 )P is equal 
to a 2 . It follows from this that p\a n if and only if (2 n — 3 )P reduces to the identity 
in E (¥ p ), and so a prime p divides a term in the Somos-4 sequence if and only if 
(0,0) £ E( F p ) has odd order. In ,8], Rafe Jones and the third author prove that 
the density of primes dividing some term of the Somos-4 sequence is + The goal 
of the present paper is to prove an analogous result for the Somos-5 sequence. 

Let 7 t' ( x ) denote the number of primes p < x so that p divides some term in the 
Somos-5 sequence. We have the following table of data. 


X 

7r'(x) 

ir'(x) 

n{x) 

10 

3 

0.750000 

10 2 

12 

0.480000 

10 3 

83 

0.494048 

10 4 

588 

0.478438 

10 5 

4539 

0.473207 

10 6 

37075 

0.472305 

10 7 

314485 

0.473209 

10 s 

2725670 

0.473087 

10 9 

24057711 

0.473134 

10 10 

215298607 

0.473129 


Our main result is the following. 


Theorem 1. We have 


lim 

x —yoo 


n'(x) 

n(x) 


5087 

10752 


0.473121. 


The Somos-5 sequence is related to the coordinates of rational points on the 
elliptic curve E : y 2 + xy = x 3 + x 2 — 2x. This curve has E{ Q) = Z x Z/2Z and 
generators are P = (2,2) (of infinite order) and Q = (0,0) (of order 2). We have 
(see Lemma [3]) that 


mP + Q = 


2 a 2 3 

^ra+2 ^m^ra+4 Q' rn CLm-\-6 ^m+2 


' t ra +2 


'*' 771+2 


It follows that a prime p divides a term in the Somos-5 sequence if and only if the 
reduction of Q modulo p is in (P) C P(F p ). Another way of stating this is the 
following: there is a 2-isogeny </>: E — > E r , where E' : y 2 + xy = x 3 + x 2 + 8x + 10 
and 

/ x 2 - 2 x 2 y + 2x + 2y\ 

= ’—?—)• 

The kernel of 4> is {0, Q}. Letting R = </>(P) we show (see Theorem [I} that p divides 
some term in the Somos-5 sequence if and only if the order of P in P(F p ) is twice 
that of R in E'(¥ p ). 

A result of Pink (see Proposition 3.2 on page 284 of m) shows that the Padic 
valuation of the order of a point P (mod p) can be determined from a suitable 
Galois representation attached to an elliptic curve. For a positive integer fc, we 
let Kk be the held obtained by adjoining to Q the x and y coordinates of all 
points Bk with 2= P. There is a Galois representation pE, 2 k '■ Gal(A\./Q) —► 
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AGL2(Z/2 fe Z) and we relate the power of 2 dividing the order of P in E{ F p ) to 
P_E, 2 fc ( <J p)) where a p is a Frobenius automorphism at p in Gal(ATfc/Q)- Using the 
isogeny <f> we are able to relate Pe 2 k ( a p) an d Pe' 2 k ~ 1 (o>), obtaining a criterion that 
indicates when p divides some term in the Somos-5 sequence. We then determine 
the image of pE, 2 k for a ll k. 

Once the image of pe 2 k is known, the problem of computing the fraction of 
elements in the image with the desired properties is quite a difficult one. We 
introduce a new and simple method for computing this fraction and apply it to 
prove Theorem |TJ 

Acknowledgements. The first and second authors thank the Wake Forest Under¬ 
graduate Research and Creative Activities Center for financial support. The authors 
used Magma JTj version 2.20-6 for computations. 


2. Background 


If E/F is an elliptic curve given in the form y 2 Faixy+a^y = x 3 -\-a 2 x 2 Pa^x+aQ, 
the set E(F) has the structure of an abelian group. Specifically, if P, Q G E(F ), 
let R = (x, y) be the third point of intersection between E and the line through P 
and Q. We define P + Q = (x, —y — a\X — 03). The multiplication by m map on 
an elliptic curve has degree m 2 , and so if E/C is an elliptic curve and a G E( C), 
then there are m 2 points j3 so that m/3 = a. 

If K/Q is a finite extension, let Ok denote the ring of algebraic integers in K. A 
prime p ramifies in K if pOk = n[=i PT and some e* > 1, where the p; are distinct 
prime ideals of Ok- 

Suppose K/Q is Galois, p is a prime number that does not ramify in K , and 
pO K = nf=i pi. For each i, there is a unique element <7 G Gal(A'/Q) for which 

cr(a) = a p (mod p*) 


for all a G Ok- This element is called the Artin symbol of p, and is denoted 


m. if Mi, 


| W/Q ’ 


and 


K/Q 

Pi 


are conjugate in Gal(AT/Q) and [—^] := 


| :!<*<g|isa conjugacy class in Gal(/\/Q). 


The key tool we will use in the proof of Theorem |T| is the Chebotarev density 
theorem. 


Theorem 2 ([7], page 143). If C C Gal(A'/Q) is a conjugacy class, then 
#{p <x:p prime , [^] = C} \C\ 


lim 

X—YOO 


r(x) 


|Gal(A7Q)| ■ 


K/Q 

P 


equally often. 


Roughly speaking, each element of Gal(/\/Q) arises as 
Let E[m\ = {P G E : mP = 0} be the set of points of order dividing m on E. 
Then Q(A[?n])/Q is Galois and Gal(Q(A[m])/Q) is isomorphic to a subgroup of 
Aut(A[m]) = GL 2 (Z/?nZ). Moreover, Proposition V.2.3 of [T3] implies that if a p is 
a Frobenius automorphism at some prime above p and r : Gal(Q(A[m])/Q) —>• 
GL 2 (Z/mZ), then tr r(a p ) = p + 1 — #E( F p ) (mod m) and det(r(cr p )) = p 
(mod m). Another useful fact is the following. If K/Q is a number field, p is 
a prime ideal in Ok above p, ged (m,p) = 1 and P G E(K)[m] is not the iden¬ 
tity, then P does not reduce to the identity in E(Ok/ p)- This is a consequence of 
Proposition VII.3.1 of IT3l . 
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or as 3 x 3 matrices 

a b 0 
c d 0 

, where v = [e f] and M = 

a b 

r r\ 


_ e / K 


L CL 


We will construct Galois representations attached to elliptic curves with images 
in AGL2(Z/2 fc Z) = (Z/2 fc Z) 2 x GL2(Z/2 fc Z). Elements of such a group can be 
thought of either as pairs (v,M), where v is a row vector, and M € GL2(Z/2 fc Z), 


In the former 


notation, the group operation is given by 

(v \, .Mi) * (v 2 , M 2 ) = (hi + 1 ’ 2 -Mi, M 2 M 1 ). 

3. Connection between the Somos-5 sequence and E 

Lemma 3. Define P = (2, 2) and Q = (0,0) on E : y 2 + xy = x 3 + x 2 — 2x. For 
all to > 0, we have the following relationship between the Somos-5 sequence and E: 

D , ^ f ^m+2 a m U m _)_4 4u m U m -)_2tt m -(-4 (-6 ^m-f-2 

mf + y I ~i - 

V u m +2 u m +2 

Proof. We will prove this by strong induction. A straightforward calculation shows 
that the base cases m = 0 and to = 1 are true. For simplicity’s sake, we will denote 
R- = a m , b — <z m _|_i, c = u m -)_ 2 , d = 3 , e = u m _)_ 4 , / = Um +5 and g = a m +e- Our 

inductive hypothesis is that 

_ _ ( c 2 — ae 4 ace — a 2 q — c 3 ' 

”- p + 0=(—■- 

We will now compute (to + 2)P + Q. 

To find the a; and y coordinates of (m + 2)P+Q, we add 2 P = (1, —1) to m.P+Q. 
If w is the slope and v is the ^-intercept, the line between 2 P and mP + Q is 
y = wx + v with w = ag ~ e 4ce and v = ~ ag ~ 3ce . Substituting this into the equation 

for E, we find the x-coordinate of 2 P + (mP + Q ) to be r x = a 9 lace ^ e ff e +9c e ■ 
A straightforward but lengthy inductive calculation shows that if 

F(a, c, e, g) = a 2 g 2 — 7 aceg + ae 3 + c 3 g + 8 c 2 e 2 , 

then P(a n , a„ + 2 , a„+ 4 , a„+ 6 ) = 0 for all n. Since F(a,c,e,g) = 0, we know that 
r x — F(a,c,e,g) = r x . Therefore, we know that r x = C9 e t e ■ 

Denote the ^-coordinate of (m + 2)P + Q as r y . We compute that r y = g( a g~ 3ce ) _ 
Using that r y = r y — , we find that r y = 4ceg ~ e 3 l ~ e . Therefore, it is evident 

that 


(to + 2 )P + Q = 


2 /i 2 3 

^ m _|_4 Q ' m + 2 ^ , m +6 4 tt rn -[_ 2 ^' m + 4 ^ , m +6 ^ m + 2^ m +8 ^ m +4 


- <, m+4 


x m+4 


□ 


Let A' be given by E' : y 2 + xy = x 3 + x 2 + 8x + 10 and let R = (1,4) G E'( 

We have a 2-isogeny f> : E —> E' given by 

,, , / x 2 - 2 x 2 y + 2x + 2y’ 

<P(x, y) = -,- 2 - 

\ X x z 

Theorem 4. If p is a prime that divides a term in the Somos-5 sequence, the order 
of P = (2,2) in E(¥ p ) is twice the order of R = (1,4) in E'(¥ p ). Otherwise, their 
order is the same. 
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Proof. If p divides a term in our sequence, say a m , we know from our previous 
lemma that the denominators (to —2)P + <3 are divisible by p. Therefore, modulo p , 
(m — 2)P+Q = 0. The point Q has order 2, so adding Q to both sides we know that 
(to — 2 )P = Q. Therefore, we can deduce that Q € (P). We have ker(0) = {<3,0} 
(see Section 3.4 of [2J). Therefore, if <f> is restricted to the subgroup generated by 
P, we have |ker(0)| = 2. Since <j>{P) = R, by the first isomorphism theorem for 
groups, TE^yj = |(-R)I- It follows that |P| = 2 ■ \R\. 

Alternatively, assume p does not divide a term in the Somos-5 sequence. So, 
there is no m such that mP + <3 = 0 modulo p, which implies that Q qL (P). 
Therefore, the kernel of f> restricted to (P) is {0} and so |P| = \4>(P)\ = |P|- □ 


4. Galois representations 


Denote by E[2 r ] the set of points on E with order dividing 2 r . Denote K r as the 
field obtained by adjoining to Q all x and y coordinates of points /3 with 2 r (3 = P. 
For a prime p that is unramified in K r , let a = 


for some prime ideal pi 
above p. Given a basis ( A,B) for P[2 r ], for any such a € Gal(iF r /Q), we have 


K r 


a (/3) = /3 + eA + fB. Also, a (A) = aA + bB and a (P) = cA + dB. Define the map 
PE, 2 k '■ Gal(/v r /Q) —► AGL 2 (Z/2 fc Z) by p E 2 it(cr) = (v,M) where M = 


and 


v = [e /]. Let r : Gal(AT r /Q) —> GL 2 (Z/2 fc Z) be given by t(ct) = M. 

Let S = {/? G E(C) : m ■ /? G E{K)} and let L be the field obtained by adjoining 
all x and y coordinates of points in S to K. Then the only primes p that ramify 
in L/K are those that divide m and those where E/K has bad reduction (see 
Proposition VIII.1.5(b) in [J3]). For E : y 2 + xy = x 3 + x 2 — 2x, the conductor of 
E is 102 = 2-3-17 and so the only primes that ramify in K r /Q are 2, 3 and 17. 

Note that, if p is unramified, there are multiple primes p, above p which could 
result in different matrices Mi and v % . However, properties we consider of these Vi 
and Mi do not depend on the specific choice of p ; . The map depends on the choice 
of basis for E[2 r ], we choose this basis as described below in Theorem [7] 

Let fi r e E(C) be a point with 2 r fi r = P. We say that fj r is an r th preimage of 
P under multiplication by 2. Let p be a prime with p / 2, 3 or 17, a = , 

and (v,M) = pE, 2 r (<x). Assume that det(/ — M) ^ 0 (mod 2 r ). This implies that 
#P(F p ) ^ 0 (mod 2 r ). 


Theorem 5. Assume the notation above. Then 2 h P has odd order in E(¥ p ) if and 
only if 2 h v is in the image of I — M. 


Proof. First, assume 2 h v is in the image of / — M. This means that x = 2 h v + xM 
for some row vector x with coordinates in (Z/2 r Z) 2 . If this is true for x = [e /], 
define C := 2 h f) r + eA + fB. We know then that cr(C) = C. We have 2 r C = 
2 h {2 r fi r ) = 2 h P. If \C\ is odd, then clearly |2 r C'| = \2 h P\ is also odd. 

If |C| is even, then every multiplication of C by 2 cuts the order by a factor of 2 
until we arrive at a point of odd order. Since |P(F p )| = det(I — M) jk 0 (mod 2 r ), 
the power of 2 dividing |Cj is also less than r, and so |2 r Cj = |2^P| is odd. 

Conversely, assume that |2 ft P| is odd. Let a be the multiplicative inverse of 2 r 
modulo \2 h P\ and define C := a2 h P G E(¥ p ). Then 2 r S = 2 h P and so 2 r (C — 
2 h [3 r ) = 0. It follows that C = 2 h /3 r + yA + zB e E(¥ p ) for some y,z € Z/2 r Z. 
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This implies that there is a Frobenius automorphism cr £ Gal(A' r /Q) for which 
cr(C) = C (mod pi) for some prime ideal p; above p. 

We claim that er(C) = C (as elements of E(K r )). Note that a(C) — C £ E[2 r ] 
and a(C) — C reduces to the identity modulo p*. Since reduction is injective on 
torsion points of order coprime to the characteristic, and p is odd, it follows that 
<j{C) = C. It follows that if pE, 2 r (u) = (v,M) then 2 h v = (/— M) [y z], which 
implies that 2 h v is in the image of / — M. □ 


The following corollary is immediate. 

Corollary 6. Let o be the smallest positive integer so that 2°v = (I — M)x for 
some x with entries in (Z/2 r Z) 2 . Then 2° is the highest power of 2 dividing |P|. 

The following theorem gives a convenient choice of basis for E[2 k ] and E'[2 k ], 

Theorem 7. Given a positive integer k, there are points A k ,B k £ E(C) that 
generate E[2 k ] and points Ck,D k £ E'(C) that generate E'{2 k ] so that 4>(A k ) = Ck 
and <f>{Bk) = 2 D k . These points also satisfy the relations: 

2Ak = Ak~i, 2Bk = Bk- 1 , 2Ck = Ck- i, and 2Dk = Dk- 1 . 

Proof. We will prove this by induction. Recall that f> : E —> E' is the isogeny with 
ker <f> = {0,T} where T = (0,0). Let </' : E' —> E be the dual isogeny, and note 
that (j) o <j>'(P) = 2P. Base Case: Let k = 1. We want to find (Ai,Bi) to generate 
E[2] and (C\,Di) to generate E'[2] so that (j){Ai ) = C/ and f>(Bi) = 2D\. We set 
B\ = (0,0), and choose A\ to be any non-identity point in E[2] other than (0,0). 
We set C\ = (j){Ai ) = (—5/4, 5/8) and choose D\ to be any non-identity point in 
E'[2] other than Ci. Note that <p'(Di) = B\. 

Inductive Hypothesis: Assume ( Ak,Bk ) = E[2 k ] and ( Ck,Dk ) = E'[2 k ] so that 
4>{A k ) = C k , <p{B k ) = 2 D k , and 4>'{D k ) = B k . Moreover, D k fL <j)(E[2 k ]). 

Since | ker^| = 2, we have that cf>(E[2 k+1 ]) A E'[2 k }. Hence, we can choose B k+ 1 
so that <j>{Bk + 1 ) = D k . Then 2B k +i = <j)'\(j>{Bk+ 1 )) = 4>'(D k ) = B k . We choose 
Dfc+i so that (j>'( y D k + 1 ) = -Bfc+i. Note that 2D k +\ = 4>(B k + 1 ) = D k and so Pfc+i £ 
E'[2 k+1 }. Now we pick A k+ \ so that 2A k +\ = A k and define Cfc+i = 4>(A k+ 1 ). 

By our Inductive Hypothesis, ( A k ,B k ) = E[2 k ], This implies that ( A k )r\(B k ) = 
0, which in turn implies that (2A k +i) fl (2B k+ i) = 0. Let C £ {A k + 1 ) (~l {B k+ 1 ). 
Then, C = aA k+ i = bB k+ i- Because \g m \ = |c| = = ^zr, 

where ord 2 (n) is the highest power of 2 dividing n, it follows that either a and b 
are both even, or they are both odd. If a and b are even, then C £ ( A k ) D ( B k ) = 0, 
which is a contradiction. If a and b are odd, then |Cj = 2 k+1 but 2C £ (Ak)r\(B k ) = 
0, which is also a contradiction. It follows that (A k + 1 ) fl {B k + 1 ) = 0, which gives 
that E[2 k+1 ] = (A k + 1 , B k +i). 

Now we show that (Ck+i, D k +\) = E'[2 k+1 ], by way of showing that (C k + i) D 
(D k +i ) = 0. We have shown that (A k +i, B k +i) = E[2 k+1 ], and so 4>{E[ 2 fe+1 ]) = 
(Cfc+i, 2D k +i). We want to show that D k+ 1 ^ 4>{E[ 2 fc+1 ]). 

If Pfc+i £ <^(P[2 fe+1 ]), then D k +\ = aCk+i-\-‘2bDk+\- So, aCk+i~\-(2b—\)D k +i = 
0. Since (26— 1) is odd, (26— l)D k +i has order dividing 2 fc+1 . Hence, aC k + i has 
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order dividing 2 fc+1 . We can then see that 

2aC k +i + 2(26 — l)Dfc_|_i = 0 
aC k + (26 - 1 )D k = 0 
=> a = (26 — 1 ) = 0 (mod 2 fc ), 

which is a contradiction. This implies that <j>(E[ 2 fc+1 ]) is an index 2 subgroup of 
(C k+ i, D k+ i) of order 2 2fc+1 , and so {C k +i , D k+ i) = E'[2 k+1 ], This proves the 
desired claim. □ 


Recall the maps : Gal(AA/Q) —> AGL 2 (Z/2 fc Z) and r : Gal(Afc/Q) —> 

GL 2 (Z/2 fc Z), defined at the beginning of this section. In |12j . an algorithm is given 
to compute the image of the 2-adic Galois representation r. Running this algorithm 
shows that the image of r (up to conjugacy) is the index 6 subgroup of GL 2 (Z/2 fc Z) 


generated by 


1 

0 


1 

1 ’ 


7 0 
2 1 


and 


5 0 
2 1 


Moreover, the subgroup generated by 


the aforementioned matrices is the unique conjugate that corresponds to the basis 
chosen in Theorem [TJ 


Theorem 8. If PE, 2 k (cr) = (v,M) where v = (e,/), then e = 0 (mod 2) if and 
only if M = 1,7 (mod 8 ). 


Proof. We will show that e = 0 (mod 2) and M = 1,7 (mod 8 ) if and only if 

<t(V2) = y/2. 

Let /3i be a point in A(A'i) so that 2/3\ = ( 2 , 2 ). We pick a basis (A±,Bi) 
according to Theorem [7] We have a(j3\) = /3 1 + eA\ + fB\, where e, / G Z/2 fc Z. 

Let <p : E —> E' be the usual isogeny and note that B i G ker </>. Thus, <f>(cr(/3i)) = 
</>(/?! + eA\ + fBi) = </>(/3i) + e(j>(Ai). It follows that e = 0 (mod 2) if and only 
if er(0(/3 1 )) = = <(>(/3i). A straightforward computation shows that the 

coordinates of cf>((3 1 ) generate Q(-\/2). It follows that e = 0 (mod 2) if and only if 
<r(y/2) = y/2. 

Finally, suppose that cr is the Artin symbol associated to a prime ideal p above 
a rational prime p. By properties of the Weil pairing (see m, Section III. 8 ), we 
have that ( 2 k = e 27 ™/ 2 * G Q(If[2 fe ]), and that a(( 2 k ) = C 2 fc t ^ M '* = C 2 fc- Since 
y/2 = Cs + Cs" 1 ! it follows easily that <r(y/2) = y/2 p = 1, 7 (mod 8 ) and hence 
cr (y/2) = y/2 if and only if det(M) = 1,7 (mod 8 ). □ 


For k > 3, define I k to be the subgroup of AGL 2 (Z/2 fc Z) whose elements are 


ordered pairs {(v,M)} where v = [e / 


group generated by 


and 


, the reduction of M mod 8 is in the 
5 01 

2 i , and e = 0 (mod 2 ) if and only 


if det(M) = 1 or 7 (mod 8 ). By Theorem [5] and the discussion preceeding it, we 
know that the image of p : Gal(AA/Q) —> AGL 2 (Z/2 fc Z) is contained in I k . 

We now aim to show that the map pEp k '■ Gal(Afc/Q) —> I k is surjective for 
k > 3. By [13; (page 105), if we have an elliptic curve E : y 2 = x 3 + Ax + B, the 
division polynomial rpm G Z [A, B, x , y\ is determined recursively by: 

-i/’i = 1, ip 2 = 2y, rp 3 = 3cc 4 + 6Ax 2 + 12Bx — A 2 , 


V >4 = 4 y(x 6 + 5 Ax 4 + 20Bx 3 - 5 A 2 x 2 - 4 ABx - 8 B 2 - A 3 ), 

1p2m+l = V’m+2'^m — ‘4’m— l^m+l i 2yip2m = V , m(V , m+ lV’m—1 — V’m—2' i /’m+l)- 
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We then define 0 m and ujm as follows: 

07n — *^0m 0771+10771—1 

42/Wm = ~ V’m-2^m+l- 

If A = — 16(4A 3 + 27B 2 ) ^ 0, then (f> m (x) and i/j m (x ) 2 are relatively prime. This 
also implies that, for P = (xo,yo) G A, 


\m\P = 


4>m(P) W ra (P) 


+ ro (P) 2 + ro (P)V ' 

Lemma 9. The map pe ,8 '■ Gal(/< 3 ,Q) —»• /3 is surjective. 


Proof. The curve E is isomorphic to E 2 '■ y 2 = x 3 — 3267x+45630. The isomorphism 
that takes E to E 2 takes P = (2,2) on E to P 2 = (87, 648) on E 2 . 

We use division polynomials to construct a polynomial fix) whose roots are 
the ^-coordinates of points 0s on E 2 so that 8 + = P 2 . By the above formulas, 
8 P 2 =(^^,|^). Since P 2 = (87,648), 

/Or) = +(P 2 )-87+(P 2 ) 2 = 0 


will yield the equation with roots that satisfy our requirement. This is a degree 64 
polynomial. By using Magma to compute the Galois Group of f(x), we find the 
order to be 8192. A simple calculation shows that J 3 has order 8192 and since f(x) 
splits in AT 3 /Q, we have that Gal(A' 3 /Q) = J 3 . □ 


To prove the surjectivity of PE, 2 k i we will consider the Frattini subgroup of Ik- 
This is the intersection of all maximal subgroups of Ik . Since Ik is a 2-group, every 
maximal subgroup is normal and has index 2. It follows from this that if g G Ik, 
then g 2 G 4>(/ fc ). 

Lemma 10. For 3 < k, <£>(/*,) contains all pairs (v, M ) such that v = 0 (mod 4) 
and M = I (mod 8 ). 


Proof. We begin by observing that for r = k, (0,1) G ^(h)- We prove the result 
by backwards induction on r. 

Inductive Hypothesis: &(h) contains all pairs (0 ,M),M = I (mod 2 r ). Let 
g = I + 2 r ~ 2 N, and let h = I + 2 r ~ 1 N. If r > 5, then a straightforward calculation 
shows that (0 ,g) G Ik■ So, (0,g ) 2 = (0 ,g 2 ) G §(Ik)- Therefore, for r > 3, 

g 2 = 1 + 2 r ~ l N + 2 2 r ~ 4 N 2 = h (mod 2 2r - 4 ). 

By the induction hypothesis, (0,g 2 fo -1 ) G $(Ik), and so (0,h) G $(/fe). 

So, for k > r > 4, all pairs (0 ,M),M = I (mod 2 r ) G $(/fc). We will now con¬ 
struct I 4 , compute $( 14 ), and determine if $( 14 ) G {(u, M) :v = 0 (mod 8 ), M = I 
(mod 8 )}. A computation with Magma shows that 


h 


'1 

1 

o' 


'7 

0 

O' 


'5 

0 

O' 

0 

1 

0 

1 

2 

1 

0 

1 

2 

1 

0 

0 

0 

1 


0 

0 

1 


1 

0 

1 


We then construct ^(h) and then </> : $( 14 ) GL 3 (Z/ 8 Z). We check that ker cj) 

has order 64 and this proves the desired claim about $( 14 ). 
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Now, observe that if v\ = (2x,2y), then {v\,1) £ Ik and so (2v±,1) = (vi ,/) 2 £ 
4>(/*,), and so &(I k ) contains all pairs (v, J) with v = 0 (mod 4). Finally, for any 
matrix M = I (mod 8 ), we have 

fa,/) * (0,Af) = ft,M) €$(/*) 

and this proves the desired claim. □ 

Finally, we prove the desired surjectivity. 

Theorem 11. XTie map PE, 2 k '■ Gal(JCfc/Q) —> Ik is surjective for all k > 3. 

Proof. Suppose to the contrary the map is not surjective. Lemma [lU] implies that 
if pE, 2 k is not surjective, the image lies in a maximal subgroup M which contains 
the kernel of the map from Ik —> 1 3 , and so the image of pe, s must lie in a maximal 
subgroup of I 3 . This contardicts Lemma [9l and shows the map is surjective. □ 

Now, we indicate the relationship between ps, 2 k an d pE\ 2 k - Let a £ Gal(ATfc/Q). 
If 0 k is chosen so 2 k j3 k = P, then 

a(A k ) = aA k + bB k , 
tx(B k ) — cA k “t - dB k , 

& (0 k) = 0k + eA k + f B k . 

Applying to these equations, we have 

<j)(a(A k )) = aC k + 2bD k = a(<j)(A k )) = a(C k ), 

<t>(<r(B k )) = cC k + 2 dD k = a(4>(B k )) = cr( 2 D k ), 

H a (0k)) = <£(&) + eC k + 2 fD k = a(<j>(0 k )) = a (0 k ), 

where 2 k /3 k = R on E'. Using the relations from Theorem 0 we have that 2D k = 
D k ~ 1 and 2 C k = C k - 1 - This gives 

&(Hk— 1 ) — aC k -i 4“ 2 bD k —i, 

Q 

<x(D k - 1 ) = —Ck—i + dD k _ 1 . 

Thus, the vector-matrix pair associated with Pe '^- 1 is (v 1 , M'), where if = [e 2/] 

and M' = . 

[2 ^. 

Let (v,M) be a vector-matrix pair in I k . Suppose that o is the smallest non¬ 
negative integer so that 2°v is in the image of (I — M). Thus there are integers C\ 
and C 2 (not necessarily unique) so that 2°v = C\X\ + C 2 X 2 , where x\ and X 2 are the 
first and second rows of / — M. 

Lemma 12. Assume that det(M — I) ^ 0 (mod 2 fc ). IfciXi+c 2 X 2 = d\X\-\-d 2 X 2 , 
then Ci = di (mod 2 ). 

Proof. The assumption on det(M — I) implies that ker (M — I) has order dividing 
2 fe_1 . However, if C\X\ + C 2 X 2 = d\X\ + d. 2 X 2 : then [ci — d\ C 2 — cfe] is an element 
of ker (M — I). If Ci ^ d\ (mod 2), then this element has order 2 fc , which is a 
contradiction. □ 

The above lemma makes it so we can speak of ci mod 2 unambigously. We now 
have the following result. 
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Theorem 13. Assume the notation above. Let o' be the smallest positive integer 
so that 2°' v' is in the image of I — M. If det (M — I) ^ 0 (mod 2 k ~ 1 ), then o ^ o' 
if and only if C\ is even. 


Proof. Let yi and j /2 be the first two rows of I — M'. A straightforward calculation 
shows that if 2 °v = C\X\ + C 2 X 2 , then 2 °v' = C\y\ + 2 c 2 y 2 - If Ci is even, then it 
follows that 2 °~ 1 if = (ci/ 2 )yi + C 21/2 and so 0 / o'. 

Conversely, if o ^ o', then o' = o — 1 and so 2°~ 1 u' = d\y\ + d 2 y 2 - We have then 
that 

2 °v = 2 dixi + d 2 X 2 (mod 2 k ~ l ). 


So if x = 


2 di 

d 2 


we have x(I — M) = 2°v (mod 2 k 1 ). If there is a vector x' with 


x ^ x’ (mod 2) so that x !(/ — M) = 2°v (mod 2 fe_1 ), then x — x' is in the kernel 
of I — M (mod 2 k ~ 1 ). However, the order of x — x! is 2 k ~ 1 and this contradicts the 
condition on the determinant. This proves the desired result. □ 


5. Proof of Theorem [T] 


Theorem [4] states that a prime p divides a term in the Somos-5 sequence if and 
only if the order of (2,2) G E( Q) is different than the order of (1,4) G E'(Q). Recall 
that o, the power of two dividing the order of P, is the smallest positive integer 
such that 2°v G im(J — M), and o' is the power of two dividing the order of R. 
Suppose that det(I — M) ^ 0 (mod 2 fc_1 ). We have 2°v G im (I — M) if and only 

^ , xi = [l — a —6], and X 2 = [—c 1 — d]. 

We know that 0 / 0 ' if and only if ci is even. For the remainder of the argument, 

a b 0 

we will consider elements of 4 as 3 x 3 matrices 


if c\X\ +C2X2 — 2°v, where M 


: 

f 


We let I - M = 


and define A = 7 / — <5e, B = af — /3e, and C = a5 — /3y. We define 


M^(fL/2 k 1,) to be the set of 3 x 3 matrices with entries in Z/2 fc Z whose third 
column is zero. We will use ord 2 (r) to denote the highest power of 2 dividing r for 
r G Z/2 fc Z. If r = 0 G Z/2 fc Z, we will interpret ord 2 (r) to have an undefined value, 
but we will declare the inequality ord 2 (r) > k to be true. 

Solving the equation C\X\ + C 2 X 2 = 2°u using Cramer’s rule gives that ciC = 
— 2°A and C 2 C = 2°B. Assuming that ci is even and o > 0 implies that C 2 must 
be odd, and this implies that ord 2 (P) < ord 2 (C). Moreover, since the power of 
2 dividing ciC must be higher than that dividing C 2 C it follows that ord 2 (P) < 
ord 2 (A). Conversely, if ord 2 (P) < ord 2 (A) and ord 2 (P) < ord 2 (C), then o > 0 and 
Ci is even. Therefore, our goal is counting of elements of Ik with ord 2 (A) > ord 2 (P) 
and ord 2 (C) > ord 2 (P). For an M 0 G M^(Z/2 r Z), define 


v(Mo,r, k) 
y(M 0 ,r) 


# {M G M 3 °( Z/2 fc Z) :M = M 0 mod(2 r ), ord 2 (A), ord 2 (C) > ord 2 (P)} , 
rj(M 0 ,r,k) 

|J 3 | -64 fe -3- 
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Roughly speaking, /x(Mo,r) is the fraction of matrices M = Mq (mod 2 r ) in Ik 
with the property that p E . 2 k (o>) = M implies that p divides a term of the Somos-5 
sequence. 

Theorem 14. We have 

lim = p(I-M, 3). 

x->oo 7 Tlx) z ' 

V ' M&I 3 

Before we start the proof, we need some lemmas. The first is straightforward, 
and we omit its proof. 

Lemma 15. If a £ Z/2 fe Z, then the number of pairs (x,y) £ (Z/2 fe Z) 2 with xy = a 
(mod 2 k ) is (ord 2 (a) + l)2 fe_1 , where if a = 0 (mod 2 k ), we take ord 2 (a) = k + 1. 

Lemma 16. The number of matrices M £ M2(Z/2 fc Z) with det(M) = 0 (mod 2 fc ) 
is 3 ■ 2 3fc " 1 - 2 2fc " 1 . 

Proof We counting quadruples ( a,b,c,d ) with aci = 5c (mod 2 fc ). By Lemma H5l 
this number is equal to 

E ((ord 2 (a) + l)2 fc_1 ) 2 , 

a£Z/2 fc Z 


which can easily be shown to equal 3 • 2 


3fc-l 


- 2 


2fe-l 


□ 


Proof of Theorem m For k > 1, let G = Gal(A^/Q) and a £ G have the property 
that a = [ E^ ] for some prime ideal p C C>K k with p D Z = (p). Assume that p is 
unramified in Kk/Q and E/¥ p has good reduction at p. Let M be the 3x3 matrix 
corresponding to p E 2 *(cr), and A , B and C be the corresponding minors of/ — M. 
Then one of three alternatives occurs: 

(a) B ^ 0 (mod 2 fe ), and a higher power of 2 divides both A and C. 

In this situation (the good case), previous results ensure that the order of P 
in E(F P ) is twice the order of R in E'(¥ p ), and hence p divides some term in the 
Somos-5 sequence. 

(b) One of A or C is not congruent to 0 mod 2 fc and the power of 2 dividing B 
is equal or higher than for A or C. 

In this situation (the bad case), previous results ensure that the order of P in 
E( F p ) is equal to the order of R in E'(¥ p ) and p does not divide any term in the 
Somos-5 sequence. 

(c) A = B = C = 0 (mod 2 fc ). 

In this situation (the inconclusive case), we do not have enough information to 
determine if p divides a term in the Somos-5 sequence or not. 

Fix e > 0 and choose a k large enough so that both of the following conditions 
are satisfied: 

W Em£1 3 V \i 3 \64 k ~ 3 ^ ~ Zmci 3 ~ 3) < e/3, and 

(ii) the fraction of elements M in 1$ with C = det(/ — M) = 0 (mod 2 fc_1 ) is 
less than e/3. (This fraction tends to zero by Lemma fTBl ) 

Let C C Ik be the collection of “good” elements of Ik and let C' be the collection 
of “good or inconclusive” elements. 

By the statements above, we have that 

E M(/-M,3)-2e/3<H 

mg / 3 I k ' 
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and 


|Jfc| Mei 3 


By the Chebotarev density theorem, we have 


//{p prime : p < x is unramified in Kk and CC} |C| 

x—toc n(x) \Ik \’ 

and the same with CL 

Let r be the number of primes that either ramify in Kk/Q or for which E/Q has 
bad reduction. Then there is a constant N so that if x > N, then 


r(x) 


r #{p prime : p < x is unramified in Kk and [ fc ^l CC} 
^ 3)- e +—-<- 

r 7rm 

Me/ 3 v ’ 

and 


4/{p prime : p < x is unramified in K/ and [—Q C'} 

ir(x) 

It follows from these inequalities that for x > N, then 

''W- £,.(/-M, 3)< e . 
mc/ 3 


< ^ p(I—M, 3)+e—- 

Mel 3 




-e < 


This proves that 


lim 


n(x) 

tt'(x) 


r(*) 


mc/ 3 


□ 


Our goal is now to compute Emc / 3 ~ M, 3). To do this, we will develop 
rules to compute p(M,r) for any matrix M £ M 3 (Z/ 2 r Z) whose third column is 
zero. Observe that p(M 0 ,r) < (m ° d 2r)} = 

Also, if all the entires in M are even, then p(M,r) = jrj/^yy, 7 ’~ 1). This allows 
us to reduce to matrices where at least one entry is odd. If M £ M^(/L/2T.) is the 
zero matrix, we have 

/1 (M,l) = i / i(M/2,0)=C E KN,1) = E V( N A)- 


64 


IVeM°( Z/2Z) 


64' 


NeM$(Z/2Z) 

N^M 


It follows that m(M, 1) = ^ EiveM 3 0 (z/ 2 Z) !)• 

NjtM 

In order to determine p(Mo,r), it is necessary to consider a matrix M £ M^,(/L/2 k T) 
and examine the behavior of matrices M' G M3(Z/2 fc+1 Z) with M' = M (mod 2 k ). 
We refer to these as 'lifts’ of M. We define A, B and C to be functions defined on 
\a /3 0] 

, given by A = 7 / — fie, B = af — fie and C = ad — /? 7 . 


a matrix M = 


7 S 0 
e / 0 
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a [3 0 

Theorem 17. Let M = 7 5 0 G Mg (Z/2 fe Z) and suppose A = B = C = 0 

_e / 0 _ 

(mod 2 fe ). 

(1) J /7 or <5 odd, then p{M, k) — 0. 

(2) If 7 and 5 are both even, but one of a, j3, e or f is odd, then p(M, k) = 

1 

6-64 fc — 1 ’ 

Proof. Suppose that M G M® (Z/2 fe Z) is a matrix with 7 or 5 odd and with A = 

C = 0 (mod 2 k ). In the case that 7 is odd, the congruences A = 0 (mod 2 k ) and 
C = 0 (mod 2 k ) imply that f = ^ (mod 2 fc ) and /3 = ^ (mod 2 fc ). We then find 

that B = af — fie = a ~ e — 0 (mod 2 fe ). It follows that none of the 
lifts of M have ord 2 (B) < min{ord 2 (A), ord 2 (C)} and so fi(M,k ) = 0. A similar 
argument applies in the case that S is odd. 

Suppose now that 7 and <5 are both even. Consider M' to be a lift of M mod 
2 k+1 . Then we have 

~a! j3' 0] \a + ai 2 k [3 + /3±2 k O' 

M' = 7 ' S' 0 = 7 + 7 i 2 fc d + <5i2 fc 0 , 

e! f oj le + ei 2 k f + f x 2 k 0 _ 

where aq, /?i, 71 , <5i, ei, /1 G F 2 . If A!, B' and C' are the values of A, B, and C 
associated to M', then 

A' = A + 2 k ( 7 l f-S ie ) (mod 2 k+1 ) 

B' = B + 2 k (aif + afi - fae - f3e 1 ) (mod 2 fc+1 ) 

C =C + 2 k (aS\ — ft'yi) (mod 2 fc+1 ). 

Suppose that e or / is odd. Then the map F® —> F| given by (or, /3±, 71 , <5i, ei, / 1 ) i-»- 
( 71 / — 5\e,otif + afi — /3ie — /3ei) is surjective. It follows that of the 64 lifts 
of M', one quarter have (A' mod 2 fc+1 , B' mod 2 fe+1 ) equal to each of (2 fc ,2 fc ), 

(0,2 fc ), (2 fc , 0) and (0,0). Moreover, if A! = 0 (mod 2 fc+1 ), then we must have 

C = 0 (mod 2 fc+1 ). This is because if e! is odd, then S' = (mod 2 fc+1 ), and 

f3' = a *~, B (mod 2 fc+1 ). Plugging these into C' = aW — (3'Y gives C = 

(mod 2 fe+1 ). Since 7 ' is even, it follows that C' = 0 (mod 2 k+l ). A similar argu¬ 
ment shows that C = 0 (mod 2 fc+1 ) if f is odd. As a consequence, of the 64 lifts of 
M , 32 have n(M ', /c+1) = 0, 16 have ord 2 (-B') < oix^A') and ord 2 (.B') < ord 2 (A'). 
For these, we have p(M', k + 1) = 9 ^ . The remainder have A! = B' = C = 0 
(mod 2 fc+1 ). It follows that 

V(M, k) = 2 _ fr^—[ ' j + XI + X )' 

M'=M (mod 2 fc+1 ) 

A'=B'=C'=ej (mod 2 fc+1 ) 

Applying the above argument repeatedly gives 

p{AI, k) = 2 _ Q ^ k l .(- + _ + ... + -'J + X M(M',fc + £). 

V 7 M'=M (mod 2 fc+i ) 

T=B'=C'=0 (mod 2 fe+f ) 
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Using the bound 0 < + £) < n , 64 it + {-i ; noting that the sum contains 16 f 

terms, and taking the limit as t —> oo yields that /z(M, k) = 2 . 64 \_i w = 

i 

6 - 64 fc -! ' 

The case when a or /3 is odd is very similar. In that case, one can show 
that the 64 lifts M' have (B 1 mod 2 fc+1 , C' mod 2 fc+1 ) divided equally between 
(2 fe , 2 k ), (0, 2 fc ), (2 fe , 0) and (0,0), and that C = 0 (mod 2 k+1 ) implies that A' = 0 
(mod 2 fc+1 ). Again, one quarter of the lifts M' have B' = 2 k (mod 2 k+1 ) and 
A' = C' = 0 (mod 2 fc+1 ), and /r(M, k) = 6 . 6 lk-i ■ □ 


Let M £ M 3 °(Z/ 8 Z) be the zero matrix. We have that /i(M, 3) = ^^(M, 1) = 
£3 ‘ ei 7 SjveM°(z/ 2 Z) MYO- Of the 63 nonzero matrices in M 3 (Z/2Z) we hnd 
that 6 have £? odd and A and C even, while 36 have A or C odd. Of the remaining 
21, there are 12 that have 7 or <5 odd, and the remaining 9 have 7 and 5 both even. 
It follows that 


63 ’ 64 2 ‘ 2 


6 + 36-0 + 12-0 + 9- - 

O 


1 1 


1 


8192 7 57344 


(Note that in the denominator of /r( N , 1) we have | 7*3 164 — 2 = 8192 • (1/4096) = 2.) 

For each of the 8191 non-identity elements M of I 3 , we divide I — M by the 
highest power of 2 dividing all of the elements, say 2 r . In 3754 cases, we have 
ord 2 (U) < ord- 2 (A) and ord 2 (l?) < ord 2 (C). For each of these, fj,(I — M, 3) = . 

In 4036 cases, we have ord 2 (l?) > ord 2 (A) or ord 2 (U) > ord 2 (C) and not all of 
A, B and C are congruent to 0 modulo 2 3 ~ r . For each of these, /x(I — M, 3) = 0. 

In 365 cases, we have A = B = C = 0 (mod 2 3_r ) and 7 and 5 are both even. 
In each of these cases, /x(7 — M, 3) = 3 ^ 8 ^ 92 by Theorem flTl 

In the remaining 36 cases, we have A = B = 0 (mod 2 3 ~ r ) and one of 7 or <5 is 
odd. By Theorem H7l / i(I — M, 3) = 0. 

It follows that 

y>(/-M, 3) = 3754.-1- — 1 1 5087 

^ v ' 8192 

mg / 3 


365- 


3•8192 57344 10752 


This concludes the proof of Theorem [l] 
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